Effective Long-Range Interactions in Confined Curved 

Dimensions 

P. Schmelcher^ 

^Zentrum fiir Optische Quantentechnologien, Universitdt Hamburg, 

Lumper Chaussee 149, 22761 Hamburg, GermanyQ 

(Dated: July 14, 2011) 

Abstract 
O, 

^N , We explore the effective long-range interaction of charged particles confined to a curved low- 

^ ' dimensional manifold using the example of a helical geometry. Opposite to the Coulomb interaction 

^^ ' in free space the confined particles experience a force which is oscillating with the distance between 

the particles. This leads to stable equilibrium configurations and correspondingly induced bound 
Q^, states whose number is tunable with the parameters of the helix. We demonstrate the existence of 
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a plethora of equilibria of few-body chains with different symmetry character that are allowed to 
freely move. An outline concerning the implications on many-body helical chains is provided. 
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I. INTRODUCTION 

Dimensionality and geometry often play a key role in physical systems and are responsible 
for many of their unusual properties. This holds in particular for condensed matter systems 



a 



where celebrated effects such as the quantum hall effect pj , high temperature superconduc- 
tivity J2| or generally the physics of quantum phase transitions jsl] are connected or even 
intrinsically tight to the underlying dimensionality. In the case of cold and ultracold quan- 
tum gases % |5| , the extensive control of the external as well as internal degrees of freedom 
of the atoms achieved during the past decades allows to systematically tune the dimension- 
ality of the confining traps thereby covering the complete crossover from zero-dimensional 
dot-like ensembles via one-dimensional gases to two-dimensional layered structures. Further- 
more, individual traps can be shaped almost arbitrarily and arrays of traps can be prepared 



and manipulated by employing magnetic and/or optical fields 



Ki. 



At the same time, the 



short-range interactions among ultracold atoms can be adjusted covering weak to strongly 
correlated systems by using magnetic or optical Feshbach resonances |9l-lll|. Long-range 
interactions of dipolar character, which can be realized by, e.g., atoms with a large magnetic 
dipole moment or heteronuclear molecules exposed to electric fields, add a plethora of novel 



phenomena both on the mean-field and strongly correlated many-body level [12J. Many 
of these effects, such as the roton-maxon spectrum and novel quantum phases in optical 
lattices, are based on the anisotropic character of the dipolar interaction. In spite of these 
intriguing findings, the underlying fundamental interactions (short-range, dipolar or even 
Coulomb interaction among e.g. ions) typically take on a very few specific appearances 
only. They are either overall repulsive or exhibit a single potential well. Asymptotically, i.e. 
for large distances, typically a power law behaviour is encountered. 

The question we pose here is whether dimensionality and geometry can cooperate in 
such a way that novel effective long-range interactions emerge. To tackle this question, we 
consider in the present work long-range (Coulomb) interacting particles that are confined 
in two out of three dimensions i.e. the particles are allowed to move on a one-dimensional 
manifold (ODM) only. The latter possesses a nontrivial geometry in the sense of a non- 
vanishing curvature with respect to its embedding into three-dimensional space. Although 
the dynamics is one-dimensional, i.e. the two confined dimensions are not accessible dy- 
namically, interactions take place via three-dimensional space. This introduces effective 




FIG. 1: Schematic illustration of the helix including the radius R of the projection onto the x,y- 
plane and the pitch h. Each pair of full and broken circles indicate the relative positions of two 
(equally charged) particles for a corresponding stable equilibrium configuration. 



one-dimensional interactions along the ODM seen by the particles which can exhibit pecu- 
liar properties, such as an oscillating short-range to long-range behaviour with an adjustable 
number of equilibrium configurations depending on the curvature and near encounters of the 
ODM. Physically speaking, one might think of a two-dimensional tightly confining trap that 
freezes the dynamics of two of the three degrees of freedom of the interacting particles and 
the remaining unconfined motional degree of freedom possesses a nonzero curvature. We 
perform a comprehensive analysis of the two-body effective interactions and dynamics and 
show the existence of few-body bound states for purely repulsive interactions in three di- 
mensions. The latter leads to a plethora of stable many-body configurations of chains in the 
ODM. In this context it is worthy to note that recently 18[ cold polar molecules confined to 
a helical optical lattice have been shown to exhibit an attractive large-distance intermolecu- 
lar interaction but the short-scale interaction is repulsive due to geometric constraints which 
leads to a zero-temperature second-order liquid-gas transition at a critical applied electric 
field. 



II. LAGRANGIAN 

We consider charged particles interacting via the Coulomb potential that are confined 
to a one-dimensional manifold. From the latter we require homogeneity and at the same 
time a nontrivial geometry in the sense of a non-vanishing curvature. One of the basic 
setups fulfilling these criteria is a helix characterized by a constant curvature and a given 
pitch (see Fig.l). The charged particles are allowed to move on the one-dimensional helix 
only but are interacting via three-dimensional space. Our starting-point is the Lagrangian 
£ = T — V using three-dimensional Cartesian coordinates r^ and corresponding velocities 
f-j where T = J2i ^'^h ^ = J2i>j 9ij/\''''i ~ ^j\ ^-nd rrii are the kinetic and potential energy 
and masses of the particles, respectively, gij are the coupling constants of the interacting 
particles. Let us now specify for each particle i its position on the helix by introducing the 
corresponding (signed) path length Sj along the helix which is connected to the Cartesian 
coordinates via Xi = Rsin (si/a), Hi = Rcos {si/a) and Zi = (/i/27r) (sj/a) where R and h are 
the radius of the circle, which is the orthogonal projection of the helix onto the (x, |/)-plane. 



and the pitch, respectively. Furthermore a = \J R^ + {h/2Tx) provides the proportionality 
factor between the path length and the azimuthal angle in the (x,|/)-plane i.e. Si = a ■ (pi. 
The Lagrangian for the one-dimensional helical motion then takes on the appearance 

N 

c=T^s^ (1) 



-E 



^ 2 

9ij 



^>J m2R^ (1 - COS (i (s, - s,))) + {^y is, - s,)' 

From the above, it is evident that the kinetic energy has retained its original Cartesian form 
and the interaction potential depends only on \si — Sj\. Therefore, besides the energy, the 
total momentum is conserved and it is possible to separate the center of mass and relative 
motion which is a consequence of the spatial homogeneity on the helix. Introducing the 
center of mass S = (X]i=i "^i-^jj /^ ^^^ (^ ~ 1) relative coordinates s, = Sj+i — Si,i = 
1, ...,N — 1 yields the following many-body Lagrangian 



'^^^s' + ^T i^Tki^-Ts. (2) 




'=' •-'= ^'(^r (E:.. i,f + 2fi2 (1 - cos (1 EL, S.)) 



where M is the total mass {gij = g and g > as well as rrii = m). In the case of two particles 
this simplifies to the integrable Lagrangian 

where s = S2 — Si and we have omitted the tilde and will, from now on, disregard the 
uniform center of mass motion. The last term of eq.([n]) establishes an 'effective' interaction 
potential Ves between the particles which are constrained to move on the helix. Note, that 
the variable s encodes the distance of the two particles on the helix and therefore replaces the 
distance of two particles in Cartesian space. The first term under the square root is remines- 
cent of the radial dependence of the Coulomb potential. The second term, however, has no 
counterpart in three-dimensional space and possesses an oscillatory behaviour. It therefore 
leads to a very unusual two-body interaction which will be analyzed in the following. 

III. TWO-BODY PROBLEM 

In spite of the purely repulsive {g > 0) Coulomb interaction in free space, Ves can, 
depending on the parameters R, h, not only be nonmonotonic with respect to its derivatives 
but even develop several local minima which qualify for a local oscillatory dynamics. To 
explore this, let us analyze the necessary conditions for the occurence of extrema of Vgs. A 
vanishing force F^ = — ^ = leads to the implicit condition 

which can geometrically be interpreted as the intersection points of a sine-function and a 
straight line. Eq.(jl]) provides the s- values which correspond to the extrema of the relative 
positions of the two equally charged particles on the helix for given parameters R, h. With 
varying parameters R, h new extrema emerge if condition (jl]) is fulfilled and the sine-function 
is tangential to the straight line. This leads to explicit conditions for the corresponding s- 
values and an implicit equation for the parameters 




(5) 



cos h H- -1 =-U-75 (6) 



In particular, the above equations require h < 2'irR, i.e. the pitch has to be smaller than 
the circumference of the projected circle in order to obtain extrema in the potential energy. 
As a result of the above conditions the number of extrema and in particular minima of the 
potential V^s can be tuned by adjusting the parameters R, h. Figure 2(a) shows the effective 
potential for different values of h. An increasing number of stable equilibrium configurations 
and consequently potential wells which lead to bound states are encountered with decreasing 
value of h. This number is in principal unbounded for /i — )■ 0. The depth of these wells 
increases with decreasing h and their position corresponding to the relative position of the 
two particles decreases slightly. Decreasing h means that we have closer near encounters of 
the helix (see Fig.l), when moving onward from a given position of the helix. It is exactly 
this spatial recurrence of the helix which leads to the occurence of the potential wells and 
which turns an overall repulsive interaction in three dimensions into an oscillatory long- 
range interaction on the helix. As a consequence two repulsively interacting particles in 3D 
possess in general a certain number of stable equilibria around which a bounded oscillatory 
motion takes place. The two particles can therefore travel in (different) bound states with 
a free center of mass motion through the helix. This is further illustrated in Fig 2(b) which 
shows the phase portrait {s,ps) belonging to the Lagrangian ([3]). Clearly visible are the 
stable fixed points, the bounded librational motion around them and the scattering winding 
trajectories coming from and going to infinity. Note that in case of attractive interactions 
the sign of Vefj has to be reversed and consequently the minimum/maximum configurations 
have to be exchanged when compared to the case g > 0. The inset in Figure 2(a) shows 
how the stable fixed points, i.e. the minima, are distributed over the helix. The relative 
phases (p = ~^~ of the positions of the two particles closely fluctuate in a regular manner 
around 1 with a slight increase of the fluctuations for increasing distance of the particles. 
This indicates that the minima and maxima of the potential energy are distributed on the 
opposite and same sides of the helical winding (see Fig. 1). The step function corresponding 
to the number of minima of Kg as a function of h is shown in the inset of Fig. 2(b). It 
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FIG. 2: Effective potential and phase space, (a) 14// as a function of ^ 



2iTa 



for R = 1 and 



g = 1. Bottom to top curves with decreasing value of h show an increasing number of stable 
equilibria and corresponding potential wells. The five lowest curves (from bottom) correspond to 
h = 5.0, 2.92, 1.89, 1.51, 1.29 and the top two curves belong to h = 0.48, 0.20. The inset shows the 
position ^^^ of the alternating sequence of minima and maxima as a function of the enumerated 
minima and maxima, (b) Phase space {s,ps = ^s) cutout {m = 2,h = 1.18). The inset shows the 
number of minima of Veff as a function of /i in a double logarithmic representation, (c) Harmonic 
frequencies of the potential wells depending on t^e minimum n (counted from the minimum closest 
to the origin). Curves from top to bottom belong to h = 0.4 up to h = 1.2 in steps of 0.1, 
respectively. 



exhibits a power law behaviour with a (negative) quadratic exponent which can be inferred 
from eq.(jl]) in the hmit h ^ 0. The emergence of new minima and maxima with decreasing 
h happens via a saddle node bifurcation cascade which, of course, obeys the same scaling 
behaviour as the number of minima. 

Figure 2(c) shows the dependence of the frequencies of the oscillatory dynamics or bound 
states around the stable equilibria on the order of the minimum (counting from the one 
closest to the origin) within an harmonic approximation for different values of the pitch 
h. The frequencies monotonically decrease for minima with increasing distance between 
the particles, i.e. the potential wells become increasingly shallower covering two orders of 
magnitude with respect to their frequencies in the given example. As one would expect, the 
frequencies are found to be overall larger for smaller values of h. The ratios of the frequencies 
belonging to different values of h increase with increasing order n of the minimum. From the 
above we can conclude that the effective interaction potential is to a large extent tunable, 
i.e. the number of minima, their frequencies and positions can be adjusted by varying the 
underlying parameters R, h. 

IV. MANY PARTICLE SYSTEMS 

Proceeding to more particle systems, the potential energy landscape Ves rapidly develops 
a large number of local minima with decreasing value of the pitch (see Fig. 3(e)). The 
two limiting cases /i — )■ and /i — )■ oo are evident. The first case corresponds to placing 
the particles in an equidistant manner on a ring and the second one to a zero curvature 
straight line and accordingly no confining wells occur, i.e. we encounter the ID Coulomb 
interaction. In between these limiting cases the interesting physics happens. Figure 3(a-d) 
shows the distribution of local minima in the $i, $2-plane with varying value of the pitch 
where $j = ^^yi = 1,2 represent the two normalized distances for the three particles. For 
h = 2.0 (Fig. 3(d)) a single equilibrium configuration of equidistantly placed particles with 
$1.2 ~ 0.68 exists. The latter means that the three particles cover 1.5 windings of the helix. 
For h = 1.5 (Fig.3(c)) already three stable equilibrium configurations and correspondingly 
bound states exist all of them being equidistant $i = $2 and covering one, two and a little 
more than three windings, respectively. For h = 1.0 (Fig. 3(b)) and h = 0.5 (Fig.3(a)) we 
encounter 14 respectively 37 equilibria in the considered coordinate regime, now many of 
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FIG. 3: (a-d) Distribution of equilibrium configurations in (<I>i, <I>2)-space (<i>j = -^) yielding 
confining potential wells and bound states for three-particles in the helix. Values for the pitch h are 
0.5 (triangles), 1.0 (stars), 1.5 (squares) and 2.0 (crosses), respectively for the different subfigures. 
(e) Three-dimensional contour plot of the potential energy for three particles for h = 1.1. 



them being of asymmetric character (see in particular Fig. 3 (a)). Fig. 3(e) illustrates the 
complexity of the potential energy landscape already manifest for a comparatively high value 
h = 1.1. Finally, table I shows all four and five particle equilibria for h = 1.5 and h = 1.1. 
For the higher value h = 1.5 only configurations symmetric with respect to the midpoint 
of the four and five-particle chain occur. For h = 1.1 however, symmetric and asymmetric 
configurations are found, the latter being of increasing or decreasing order or even oscillating 
with respect to the values of $j, i = 1,2, 3. 

Extrapolating from a few to many-particle systems the potential landscape of ions in a 
helical chain will exhibit a plethora of local minima reminescent of spin glasses [13| which are 
characterized by frustrated interactions and stochastic disorder. The metastable structures 
in spin glasses are replaced here by the metastable ionic pair configurations where two ions 
are located at different windings of the helix and separated by a finite potential barrier. With 
increasing distance (in terms of windings) between the particles this barrier lowers and the 
many-particle quantum helical chain therefore comprises many metastable structures with 
largely varying time-scales. These facts will presumably lead to intriguing thermodynamical 
properties of our helical chain of ions. The non-ergodic dynamics of spin glasses for low 
temperature occurs due to trapping in the deep tales of the hierarchically frustrated energy 
landscape. In our case of the ionic helical chain we also expect non ergodic behaviour. 
However, in contrast to spin glasses the minima and corresponding metastable states follow 
a hierarchical scheme. Indeed, Monte-Carlo simulations of the statistical finite-temperature 
properties of the classical chain indicate already a very unusual low-temperature behaviour 
below and at the transition point where the ions can pass the corresponding potential barriers 
i.e. traverse the blocked windings of the chain. However, a detailed study of these properties 
goes beyond the scope of the present work and will be presented elsewhere IJ] . 

V. CONCLUSIONS 

By using the example of the Coulomb interaction in three dimensions we have shown that 
a confinement of the motion to a one- dimensional curved manifold introduces a novel type of 
effective interaction among the particles. The latter possesses a peculiar form unknown from 
fundamental interactions in Cartesian space and is to a large extent tunable by changing the 
curvature or pitch of the underlying helical configuration. Due to the helical winding, the 
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repulsive Coulomb interaction can develop for sufficiently small pitch an arbitrary number of 
local equilibrium configurations and corresponding bound states which have been explored 
in detail for few-body systems. It is evident that for many-particle chains there exist a 
plethora of local equilibria covering a wide regime of interparticle separations. It is then 
a natural step to develop the quantum physics of the many-body interacting helical chain. 
This will lead to novel structural properties such as enriched phase diagrams as well as an 
intriguing dynamical behaviour of such a chain. Open questions include the metastability 
of the observed equilibria and the quest for the existence of long-range versus short-range 
order in these systems. The analogy to spin glasses with their frustrated interactions and 
many metastable states on varying time scales suggests an equally intricate statistical finite 
temperature behaviour to occur for our many-body helical chain. 

Finally, let us briefly address the question of possible experimental realizations of the 



helix topology. Ion trap technology and in particular microfabrication of ion traps [15| offer a 
wide range of possibilities to control even individual ions and manipulate them with electric 
fields as well as optical and microwave radiation. Still, the possibility to create a helical 
configuration in this framework represents an experimental challenge. It can however be 
realized in the context of free-standing nano-objects based on few monolayers thick scrolled 
heterostructures [16| leading to helices where both the radius and pitch are of the order of 
10-100 nm. For neutral species, on the other hand side, optical interference based chiral 



microstructures can be fabricated using holographic lithography [17|. Also the work in 



ref . [18|] clearly demonstrates the experimental relevance and realizability as well as the novel 
physics of dipolarly interacting molecules in helical optical lattices. An intriguing perspective 
is here the use of evanescent fields surrounding optical nanofibers which allow to trap and 



coherently process neutral atoms or polar molecules 19 1. 
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N = A 

$1 - $2 



$, 



N = 5 

$1 - $2 



- $0 - <J). 



/l= 1.1 



0.273 
0.303 
0.471 
0.750 
0.681 
0.667 
1.766 
1.754 
1.635 



0.267 
0.473 
0.292 
0.473 
0.733 
0.689 
0.689 
0.632 
1.727 



0.273 
0.750 
0.471 
0.303 
0.681 
1.766 
0.667 
1.754 
1.635 



0.261 
0.274 
0.361 
0.328 
0.418 
0.483 
0.557 
0.669 
0.737 
0.780 
0.637 
0.307 
0.681 
1.756 
1.644 
1.790 
1.757 



0.357 
0.574 
0.239 
0.458 
0.392 
0.268 
0.256 
0.754 
0.553 
0.416 
0.790 
0.477 
0.743 
0.715 
1.739 
0.699 
0.768 



0.256 
0.574 
0.239 
0.790 
0.392 
0.553 
0.357 
0.754 
0.268 
0.417 
0.458 
0.768 
0.715 
0.743 
1.739 
0.699 
0.477 



0.557 
0.273 
0.361 
0.637 
0.418 
0.737 
0.261 
0.669 
0.483 
0.780 
0.328 
1.757 
1.756 
0.681 
1.644 
1.790 
0.307 



h= 1.5 



0.470 - 0.327 
0.672 - 0.748 



0.470 
0.672 



0.670 - 0.768 - 0.768 - 0.670 



TABLE I: Four and five-particle equilibrium configura- 
tions corresponding to bound states on the helix for two 
different values of the pitch. Positions on the helix are 
provided in terms of $« = -^^ 
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